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Abstract 

The HyperCP collaboration has recently reported the observation of three events for the decay 
S + — > pu + u~. They have suggested that new physics may be required to understand the implied de- 
cay rate and the observed distribution. Motivated by this result, we re-examine this mode within the 
standard model, considering both the short-distance and long-distance contributions. The long-distance 
part depends on four complex form-factors. We determine their imaginary parts from unitarity, fix two of 
the real parts from the S + — ► y>7 measurements, and estimate the other two with vector-meson-dominance 
models. Taking into account constraints from S + — ► pe + e~, we find that S + — > pu + u~ is long-distance 
dominated and its rate falls within the range suggested by the HyperCP measurement. 
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I. INTRODUCTION 



Three events for the decay mode S + — > p/j, + fi~ have been recently observed by the HyperCP 
(E871) collaboration [1] with results that suggest new physics may be needed to explain them. In 
this paper we re-examine this mode [2] within the standard model. 

There are short- and long-distance contributions to this decay. In the standard model (SM), 
the leading short- distance contribution comes from the Z-penguin and box diagrams, as well as 
the electromagnetic penguin with the photon connected to the dimuon pair [3]. We find that this 
contribution yields a branching ratio of order 10~ 12 , which is much smaller than the central ex- 
perimental value of 8.6 x 10~ 8 reported by HyperCP [1]. It is well known that the long-distance 
contribution to the weak radiative mode S + — > jrf is much larger than the short- distance contri- 
bution. It is therefore also possible to have enhanced long-distance contributions to S + — > pfi + fi~ 
via an intermediate virtual photon from S + — > jry. We find that the resulting branching ratio is 
in agreement with the measured value. There is, of course, still the possibility [4] that new physics 
is responsible for the observed branching ratio of S + — > py and hence that of S + — > pfi + fi~ . This 
implies that it is essential to have an up-to-date estimate of the standard-model contributions, on 
which we concentrate in this work. 

In Sec. II we update the estimate of the short-distance amplitude. We use the standard effective 
Hamiltonian for the s — > d£ + £~ transition [3] supplemented with hadronic matrix elements for 
the relevant currents. In Sec. Ill we study the long-distance contributions mediated by a real or a 
virtual photon. These can be parameterized by four (complex) gauge- invariant form-factors [2] . We 
determine the imaginary parts of these form factors from unitarity. The real parts of two of the form 
factors can be reasonably assumed to be constant as a first approximation and can then be extracted 
from the measured rate and asymmetry parameter for S + — > jry up to a fourfold ambiguity. The 
real parts of the two remaining form-factors cannot be extracted from experiment at present, and so 
we estimate them using vector-meson-dominance models. Finally, in Sec. IV we combine all these 
results to present the predictions for the rates and spectra of the two modes S + — > pfi + fi~ , pe + e~ . 
Before concluding, we discuss the implications of our analysis for the possibility that new physics 
could be present in the recent measurement by HyperCP. 



II. SHORT-DISTANCE CONTRIBUTIONS 



The short-distance effective Hamiltonian responsible for S + — > p£ + £ contains contributions 
originating from the Z-penguin, box, and electromagnetic-penguin diagrams. It is given by [3, 5] 

Ks = %OL [{*rv + ry 7 v)Orv + ry 7A 7A ] + ^ E ' W 
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where V kl are the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [6], z, y, and c are 
the Wilson coefficients, r = -V* d Vj (V* d V us ) , and 

7V = drf{\ - l5 )st % £ + , 7A = drf(\ - j 5 )s£- % j 5 £ + , 

7l = da^F^ [m s (l + j 5 ) + m d (l - 7b )] s , (2) 

with F^ v being the photon field-strength tensor. The contribution of 7l to S + — > p£ + £~ occurs via 
the photon converting to a lepton pair. The total short- distance contribution to the E + — > p£ + £~ 
amplitude is then given by 



M^+^ptr) = (p£ + r\H cH \z + ) (3) 

G p 

71 



% { V* d V us [(z 7V + ry 7V )(p\d^(l - 75 ) S |£+>r 7 / + + Ty 7 M*f(l ~ 7 5 )s|£ + K~7^ + ] 



- E Wis K™. + m d ){p\da^q v s\E + ) + (m s - m d )(p\da^q vl ,s\^)} tj^ 
j 

where q = p E — p p . 

To obtain the corresponding branching ratio, one needs to know the hadronic matrix elements. 
Employing the leading-order strong Lagrangian in chiral perturbation theory (xPT), given in 
Eq. (Bl), we find 

(p\d^s\J: + ) = -p 7 M S, (p\dn 5 s\Z + ) = {D-F)pry^ 6 E, (4) 

where D = 0.80 and F = 0.46 from fitting to hyperon semileptonic decays, and using quark-model 
results [7] we obtain 

(p\da^s\L + ) = c a pa^ , (p\da^ l5 s\Z + ) = c a pa^ 5 Z , (5) 

where c a = —1/3. Furthermore, we adopt the CKM-matrix elements given in Ref. [8], the typical 
Wilson coefficients obtained in the literature [3, 5], namely z 7V = —0.046a, y 7V = 0.735a, y 7A = 
—0.700a [3], and cL being dominated by c 7 = 0.13 [5], and the quark masses m d = 9MeV and 
m s = 120 MeV. 

The resulting branching ratio for S + — > pfi + fi~ is about 10 -12 , which is way below the observed 
value. There are uncertainties in the hadronic matrix elements, the Wilson coefficients, and the 
CKM-matrix elements, but these uncertainties will not change this result by orders of magnitude. 
We therefore conclude that in the SM the short-distance contribution is too small to explain the 
HyperCP data on S + — > pfi + fi~ . 

Now, a large branching ratio for S + — > p£ + £~ may be related to the large observed branching 
ratio for S + — > jry, compared with their respective short-distance contributions. With only the 
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short-distance contribution to S + — > pry within the SM, the branching ratio is predicted to be 
much smaller than the experimental value [4]. However, beyond the SM it is possible to have 
an enhanced short-distance contribution to S + — > jry [4] which would enhance the amplitude for 
E + — > pfi + fi~. The origin of the enhancement may be from new interactions such as W^-Wr 
mixing in left-right symmetric models and left-right squark mixing in supersymmetric models [4]. 
These types of interactions have small effects on other related flavor- changing processes such as 
K°-K° mixing, but can have large effects on S + — > jrf and therefore also on S + — > p£ + £~ . Thus 
the observed branching ratio for S + — > pj can be reproduced even if one assumes that there is 
only the short-distance contribution. More likely, however, the enhancement is due to long-distance 
contributions within the SM. In the next section we present the most complete estimate possible at 
present for these long-distance contributions. 



III. LONG-DISTANCE CONTRIBUTIONS 

In this section we deal with the contributions to S + — > p£ + £~ that are mediated by a photon. 
For a real intermediate photon there are two form factors that can be extracted from the weak 
radiative hyperon decay B :i — > B^ and are usually parameterized by the effective Lagrangian 

£ = ^ B f (a + 675) o^B, . (6) 

The two form factors, a and b, are related to the width and decay distribution of the radiative decay 
by 

Ti^Bji) = ^(|a| 2 +|6|> 3 , (7) 

7T 

dT n 2 Re 

1 + a cos , a 



dcos6 ' |a| 2 + \b\ 2 ' 

where u is the photon energy, and 9 is the angle between the spin of B { and the three-momentum 
of Bj. The measured values for S + — > are [8] 

r(£ + ^p 7 ) = (10.1 ±0.4) x KT 15 MeV , a = -0.76 ± 0.08 . (9) 

When the photon is a virtual one, there are two additional form-factors, and the total amplitude 
can be parameterized as 

M(B, t -> Bfj*) = -eG F B f [1^(0 + 675) + (?V-?7)(c + ^)] ^4 , (10) 
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where q is the photon four-momentum. We note that the a and c (b and d) terms are parity 
conserving (violating). The corresponding amplitude for B i — > B^£ + £~ is then 

-ie 2 G 

M(Bi -> BfPT) = B f (a + 6 75 ) a^B t 

- e 2 G F B fllM (c + d l5 )B i r^£+ , (11) 

where now q = p e+ +p t ~. In general a, b, c, and d depend on q 2 , and for S + — > pj* the first two 
are constrained at q 2 = by the data in Eq. (9) as 

|a(0)| 2 + |6(0)| 2 = (15.0 ±0.3) 2 MeV 2 , 

Re (a(0) 6*(0)) = (—85.3 ± 9.6) MeV 2 . (12) 
These form factors are related to the ones in Ref. [2] by 

a = lib, , b = 2ib 2 , c = ^± , d = . (13) 

As we will estimate later on, these form factors have fairly mild g 2 -dependence. If they are taken 
to be constant, by integrating numerically over phase space we can determine the branching ratios 
of S + — > p£ + £~ to be, with a and b in MeV, 

£(£+ -> p/i + /0 = [2.00 (H 2 + |6| 2 ) - 1.60 (H 2 - |6| 2 )] x lO" 10 
+ (1.05 |c| 2 + 18.2 \d\ 2 ) x 10~ 6 

+ [0.29 Re(ac*) - 16.1 Re (bd*)) x 10~ 8 , (14a) 

S(S + -> pe+e") = [4.22 (|a| 2 + |6| 2 ) - 0.21 (|a| 2 - |6| 2 )] x 10" 8 
+ (5.38 |c| 2 + 15.9 \d\ 2 ) x 10~ 5 

+ [1.51 Re (ac*) - 21.1 Re {bd*)) x 10" 7 . (14b) 

If the form factors have g 2 -dependence, the expression is different, and the rate should be calculated 
with the formula which we give in Appendix A. 

A. Imaginary parts of the form factors from unitarity 

The form factors which contribute to the weak radiative hyperon decays have been studied 
in chiral perturbation theory [9-11]. The imaginary parts of a and b for S + — > have been 
determined from unitarity with different results in the literature. Neufeld [9] employed relativistic 
baryon xPT to find, for q 2 = 0, 

lmo(0) = 2.60 MeV , Imfe(O) = -1.46 MeV (15) 
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in the notation of Eq. (6), whereas Jenkins et al. [10] using the heavy-baryon formulation obtained 

Ima(O) = 6.18 MeV , Im6(0) = -0.53 MeV . (16) 

Because of this disagreement, and since we also need the imaginary parts of the form factors c 
and d, we repeat here the unitarity calculation employing both the relativistic and heavy baryon 
approaches. 

Our strategy to derive the imaginary parts of the four form-factors in Eq. (11) from unitarity 
is illustrated in Fig. 1. As the figure shows, these imaginary parts can be determined from the 
amplitudes for the weak nonleptonic decays S + — > pir° and S + — > nn + (the vertex indicated by 
a square in Fig. 1) as well as the reactions Ntt — > N'-f* (the vertex indicated by a blob in Fig. 1). 
The weak decays have been measured [8], and we express their amplitudes as 1 

M(Y> + Ntt) = %G F m\ + N (A Nw - B Nwl5 ) £ , (17) 

where 

A nn + = 0.06 , B n7T+ = 18.53 , 

A p7T o = -1.43 , B p7T o = 11.74 . (18) 
Following Refs. [9, 10], we adopt the Ntt — > p-y* amplitudes derived in lowest-order x?T. 




FIG. 1: Unitarity cut. 



We present the details of our unitarity calculation in Appendix B. The results in the relativistic 
and heavy baryon approaches are given in Eqs. (B2) and (B8), respectively. In Fig. 2 we display 
the two sets of form factors for < q 2 < (m 2 — m N ) 2 . We note that, although only the S + — > 
nn + transition contributes to the heavy-baryon form-factors at leading order, the sizable difference 
between the Im a, or Im c, curves arises mainly from relativistic corrections, which reduce the heavy- 
baryon numbers by about 50%. On the other hand, the difference between the Im6, or Im<i, curves 
is due not only to relativistic corrections, but also to A n7T + being much smaller than A p7T o. 



We have taken the nonzero elements of 7 5 to be positive. 



6 




FIG. 2: Imaginary parts of the form factors in S + — ► p"f* , obtained using heavy baryon xPT (solid lines) 
and relativistic baryon xPT (dashed lines). 



To compare with the numbers in Eqs. (15) and (16) calculated in earlier work, we find from the 
relativistic formulas in Eq. (B2) 

Ima(O) = 2.84 MeV , Im6(0) = -1.83 MeV , (19) 

and from the heavy-baryon results in Eq. (B8) 

Ima(O) = 6.84 MeV, lm&(0) = -0.54 MeV . (20) 

Thus our relativistic results are close to those in Eq. (15), from Ref. [9], and our heavy-baryon 
numbers to those in Eq. (16), from Ref. [10]. 2 These two sets of numbers are different for the 
reasons mentioned in the preceding paragraph. 



2 Our heavy-baryon expressions for Ima(O) and Im6(0) are identical to those in Ref. [10], except that their Imo(O) 
formula has one of the overall factors of l/(m s — Tnjy) apparently coming from their approximating [(m s — 
m N ) 2 — m 2 ] 1 / 2 as to s — m N . This is the main reason for the value of Ima(O) in Eq. (16) being smaller than that 
in Eq. (20). 
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B. Real parts of the form factors 



The real parts of the form factors cannot be completely predicted at present from experimental 
input alone. For Rea(g 2 ) and Re6(g 2 ), the values at q 2 = can be extracted from Eq. (12) after 
using Eq. (19) or (20) for the imaginary parts. Thus the relativistic numbers in Eq. (19) lead to 
the four sets of solutions 

Rea(0) = ±13.3 MeV , Re 6(0) = ^6.0 MeV , 
Rea(0) = ±6.0 MeV , Re 6(0) = =f13.3 MeV , (21) 

while the heavy-baryon results in Eq. (20) imply 

Rea(0) = ±11.1 MeV , Re 6(0) = ^7.3 MeV , 
Rea(0) = ±7.3 MeV , Re6(0) = =f11.1 MeV . (22) 

Since these numbers still cannot be predicted reliably within the framework of %PT [9, 10], we will 
assume that 

Rea(g 2 ) = Rea(0) , Re6(g 2 ) = Re 6(0) , (23) 

where the q 2 = values are those in Eqs. (21) and (22) in the respective approaches. This 
assumption is also reasonable in view of the fairly mild g 2 -dependence of the imaginary parts seen 
in Fig. 2, and of the real parts of c and d below. In predicting the S + — > p£ + £~ rates in the 
following section, we will use the 8 sets of possible solutions in Eqs. (21) and (22). 

The real parts of c and d cannot be extracted from experiment at present. Our interest here, 
however, is in predicting the SM contribution, and therefore we need to estimate them. To do so, we 
employ a vector-meson-dominance assumption, presenting the details in Appendix C. The results 
for Rec(g 2 ) and Re<i(g 2 ) are given in Eqs. (C3) and (C5), respectively. In Fig. 3 we display the two 
form factors for < q 2 < (m E — m N ) 2 . We can see from Figs. 2 and 3 that c is dominated by its 
imaginary part, but that d is mostly real. 



IV. RESULTS AND CONCLUSIONS 



We can now evaluate the rates and spectra of S + — > p£ + £r resulting from the various standard- 
model contributions. Since the short-distance contributions discussed in Sec. II are very small, we 
shall neglect them. Consequently, the rates are determined by the various form factors in S + — > p^y* 
calculated in the preceding section and applied in Eq. (Al). 

In Table I, we have collected the branching ratios of S + — > pfi + fi~ and S + — > pe + e~ corre- 
sponding to the 8 sets of solutions in Eqs. (21) and (22), under the assumption of Eq. (23) for Re a 
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q 2 (GcV 2 ) q 2 (GcV 2 ) 



FIG. 3: Real parts of c and d. 

and Refc. The real parts of c and d in Eqs. (C3) and (C5) are used in all the unbracketed branching 
ratios. For the imaginary parts of the form factors, the expressions in Eq. (B2) [Eq. (B8)] contribute 
to the unbracketed branching ratios in the upper (lower) half of this table. Within each pair of 
square brackets, the first number is the branching ratio obtained without contributions from both 
c and d, whereas the second number is the branching ratio calculated with only the real parts of all 
the form factors. 

TABLE I: Branching ratios of S + — ► p(i + fj>~,pe + e~ in the standard model. The unbracketed branching 
ratios receive contributions from all the form factors, with the expressions in Eq. (B2) [Eq. (B8)] for the 
imaginary parts contributing to the numbers in the first (last) four rows. Within each pair of square 
brackets, the first number has been obtained with c = d = 0, and the second with only the real parts of 
all the form factors. 



Re a (MeV) 


Re b (MeV) 


10 8 £(£ 


f ) 


10 6 B(£- 


^ — > pe + e ) 


13.3 


-6.0 


1.6 


[2.2, 1.3] 


9.1 


[9.2, 8.6] 


-13.3 


6.0 


3.4 


[2.2, 3.1] 


9.4 


[9.2, 8.8] 


6.0 


-13.3 


5.1 


[6.7, 4.7] 


9.6 


[9.8, 9.0] 


-6.0 


13.3 


9.0 


[6.7, 8.6] 


10.1 


[9.8, 9.5] 


11.1 


-7.3 


2.3 


[2.9, 1.5] 


9.3 


[9.3, 7.2] 


-11.1 


7.3 


4.5 


[2.9, 3.7] 


9.6 


[9.3, 7.5] 


7.3 


-11.1 


4.0 


[5.1, 3.2] 


9.5 


[9.6, 7.4] 


-7.3 


11.1 


7.3 


[5.1, 6.4] 


10.0 


[9.6, 7.8] 



In Fig. 4 we show the invariant- mass distributions of the fi + fi~ pair, with M = yg 2 , that 
correspond to the smallest and largest rates of S + — > pfi + fi~ listed in Table I for both the relativistic 
baryon [(a) and (b)] and heavy baryon [(c) and (d)] cases. For S + — > pe + e~ , the mass distributions 
of the e + e~ pair, two of which are displayed in Fig. 5, differ very little from each other and are 
strongly peaked at low M ee = \f(f. Also shown in the figures are the distributions obtained with 
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Re a = 13.3 MeV 
Rc& = -6.0 MeV 



> 



+ 
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220 240 260 
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FIG. 4: Invariant-mass distributions of the lepton pair in E + — ► pjjL + \x~ corresponding to the smallest 
and largest branching ratios for the (a,b) relativistic and (c,d) heavy baryon cases in Table I. In all 
distribution figures, each solid curve receives contributions from all the form factors, each dashed curve 
has been obtained with c = d = 0, and each dotted curve involves no imaginary parts of all the form 
factors. 



c = d = (dashed curves), as well as those without contributions from the imaginary parts of all 
the form factors (dotted curves). 

We can see from Table I, Fig. 4, and Fig. 5 that the effect of the c and d contributions on the 
total rates can be up to nearly 40% in S + — > pfi + fi~, but it is much smaller in S + — > pe + e~. 
Furthermore, the contributions of the imaginary parts of the form factors can be as large as 35% 
to the p/j, + fj,~ rate and roughly 20% to the pe + e~ rate. This implies that a careful analysis of 
experimental results, especially in the case of S + — > pfi + fi~ , should take into account the imaginary 
parts of the form factors. 

For S + — > pn + /j~, HyperCP measured the branching ratio to be (8.6I54 ± 5.5) x 10 -8 [1]. It is 
evident that all the predictions in Table I for the pfi + fi~ mode corresponding to the different sets 
of form factors fall within the experimental range. For S + — ► pe + e~ , the branching ratio can be 
inferred from the experimental results given in Ref. [12], which reported the width ratio T(S + — > 
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FIG. 5: Low- mass portion of the invariant-mass distributions of the lepton pair in S + — ► pe + e corre- 
sponding to two of the branching ratios in Table I, for the (a) relativistic and (b) heavy baryon cases. 



pe + e~)/r(S + — > pir ) = (1.5 ± 0.9) x 10~ 5 and interpreted the observed events as proceeding from 
E + — > p7*, based on the very low invariant-masses of the e + e~ pair. 3 This number, in conjunction 
with the current data on S + — > pn [8], translates into i3(E + — > pe + e~) = (7.7 ± 4.6) x 10~ 6 . 
Clearly, the results for the pe + e~ mode in Table I are well within the experimentally allowed 
range. Based on the numbers in Table I, we may then conclude that within the standard model 

1.6 x lCT 8 < < 9.0 x lCT 8 , 

(24) 

9.1 x 10~ 6 < -»• pe + e~) < 10.1 x 10~ 6 . 

The agreement above between the predicted and observed rates of S + — > p£ + £~ indicates 
that these decays are dominated by long-distance contributions. However, the predicted range 
for S(S + — > pfi + fi~) is sufficiently wide that we cannot rule out the possibility of a new-physics 
contribution of the type suggested by HyperCP [1] . Motivated by the narrow distribution of dimuon 
masses of the events they observed, they proposed that the decay could proceed via a new inter- 
mediate particle of mass ~ 214 MeV, with a branching ratio of (3.ltH ± 1.5) x 10~ 8 [1]. For this 
hypothesis to be realized, however, the new physics would have to dominate the decay. It will be 
interesting to see if this hypothesis will be confirmed by future measurements. 

Finally, we observe that the smaller numbers i3(S + — > p[i + fi~) ~ 2 x 10~ 8 in Table I correspond 
to the mass distributions peaking at lower masses, ~ 220 MeV, in Fig. 4. It is perhaps 

not coincidental that these numbers are similar to the branching ratio and new-particle mass, 
respectively, in the HyperCP hypothesis above. This may be another indication that it is not 
necessary to invoke new physics to explain the HyperCP results. 



3 We note that the upper limit of 7 x 10 6 quoted in Ref. [8] and obtained in Rcf. [12] is for the presence of weak 
neutral currents in S + — > pe + e~ and not for the branching ratio of this mode. 
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APPENDIX A: DIFFERENTIAL RATE OF S" 



p£+t 



If the form factors have g 2 -dependence, before integrating over phase space to obtain the branch- 
ing ratio we should use 
c/r(£+ -> p£+£-) 



dq 2 dt 



a 2 G 2 F 
47rm| 



x I [(2m 2 + q 2 )((m p - m E ) 2 - g 2 )(m E + m p ) 2 + 2q 2 f(m p ,m^,m h q 2 ,t)]- 



t 

+ [(2m 2 + g 2 )((m p + m E ) 2 - g 2 )(m E - m p ) 2 + 2g 2 f(m p , m E , m h q 2 , t)] 

+ [(2m 2 + g 2 )((m p -m s ) 2 -g 2 ) - 2f(m p , m E , m u g 2 , t)] \c\ 2 
+ [(2m 2 + q 2 )({m p + m^) 2 - q 2 ) - 2f{m pJ m J ,,m l ,q 2 ,t)} \d\ 2 

+ 2(m E + m p )(2m 2 + q 2 ) [(m p - m E ) 2 - q 2 )] 

, \/r\2 2\ l~/ \2 21 (bd 

- 2(m E - m p )(2m l + q ) [(m p + m E ) - g J — 

where t = (p E — p e -) 2 and 

/(m p , m E , m ; , g 2 , t) = mf + (m 2 + m| - g 2 - 2t)m 2 + m 2 m| - (m 2 + m|)t + (g 2 + t)t 
with the integration intervals given by 



(Al) 



max.mm 



ml + ml + 2m 2 - g 2 ± J 1 



4m? 



m| — m 2 — g 2 ) 2 — 4m 2 g 2 



(A2) 



?L = ( m s- m P ) 2 



It is worth mentioning that, since the form factors belong to the S + — > p-y* amplitude, they do not 
depend on t. 



APPENDIX B: IMAGINARY PARTS OF FORM FACTORS IN %PT 

The chiral Lagrangian for the interactions of the lowest-lying mesons and baryons is written down 
in terms of the lightest meson-octet and baryon-octet fields, which are collected into 3x3 matrices 
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FIG. 6: Leading-order diagrams for Nir — > jry* reactions. 



(f and B, respectively [13]. The mesons enter through the exponential £ = £ 2 = exp(itp/f), where 
f — f n — 92.4 MeV is the pion decay constant. In the relativistic baryon ;\TT, the lowest-order 
strong Lagrangian is given by [13] 

4 = (Btr{d^B+ [V^B])) + m (BB)+D(B 1 ^ 5 {A^B}} + F(Bj^ 5 [A IM B]) , (Bl) 

where (•••)= Tr(- •• ) in flavor space, m is the baryon mass in the chiral limit, V" = ~(£d M £ f + 
e f + I e A" + CQ^) , and A* = f (£ 3"^ - e f <9 M £) + | eA" (f tQf - e<5e f ) , with A" being 
the photon field and Q = diag(2, — 1, — 1)/3 the quark-charge matrix. 4 The parameters D and F 
will enter our results below only through the combination D + F = 1.26. 

From C s we derive two sets of diagrams, shown in Fig. 6, which represent the Ntt p^f* 
reactions involved in the unitarity calculation of the imaginary parts of the form factors a, b, c, and 
d. It then follows from Fig. 1 that the first set of diagrams is associated with the weak transition 
S + — > nn + , and the second with S + — > pn°. Consequently, we express our results as 

where JF + (jF Q ) comes from the nir + (pn°) contribution, and write them in terms of the weak 
amplitudes A + = A n7T +, A Q = A p7T o, B + = B n7T +, and B = B pn o given in Eq. (18). Working in 
the S + rest-frame, which implies that the energies and momenta of the photon and proton in the 
final state and of the pion in the intermediate are fixed by kinematics, we define 

(B3) 



2^£ 7 - 2 


\Pn\ 


\Pj 


\-Q 2 \ 


2£.£ 7 + 2 


Pj 


\Py 


-q 2 ) 





-2|pJ 


\P P \ 


2 

— mt 

TV 


^E p 


+ 2|pJ 


\P P \ 


— ml 



4 Under a chiral transformation, B -> UBU\ B -> UBU\ -> UV*W + id^UW, and A* -> UA^W, where 
U is defined by f -> Lf £/t = f/£Rt. 
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The expression for T from each set of diagrams can then be written as 





2|pJ 


Ip 7 I/S + M«+,o)^ 




2ml 


Pi 


[(m s - m N ) 2 - q 2 \ [(m E + m N ) 2 - q 2 


2 



_ ^+,c 


m 7V 


2|P.| 




, + ln(z +>0 ) g% 




2m| 


P 7 


[(m s - mj 2 - g 2 


2 [(m s + m ff ) 2 - q 2 } 





5 + 


,0 m Af 




2|pJ|p 7 |/$ + ]ii(*_ 


1 o (c) ' 
-,oJ y+,o 




2m| 


|P 7 I 
-A 


{m N - m s ) 
+,o m w 




-m^-g 2 ] [(m s + 
2|pJ|p 7 |/!V ln (^ 


m N ) 2 - 

) a {d) ' 
-,o) y+,o 


-g 2 ] 2 


2m| 


|P 7 I 


(m s + m^) 


[(™E 


- m N ) 2 - q 2 } 2 [(m s - 


h m^) 2 


-q 2 ] 



where 

/| a) = m N m^ + (g 2 + 2m 2 + m 2 N ) - m N (3q 2 - 3m 2 + 2m 2 N ) m| 

- (g 4 - 5m 2 g 2 + 2771^ + (g 2 + m 2 ) m^) m| 

+ m N (m 2 N - q 2 ) (2g 2 - 3m 2 + m 2 N ) m s + (g 2 - m 2 N ) 2 (m 2 N - m 2 ) , 
g+ = (m N q 6 + (m N (2m N - m E ) (m N + m E ) - m 2 (377% + m s )) g 4 

+ m 2 (3m 2 - 4m^) (m w + m s ) g 2 + m 2 (m^ - m E ) 2 (m w + m s ) 3 ) , 
fo a) = 3m Ar m| - (g 2 - 2m 2 - 3m 2 N ) m| - m N (4m 2 N - 3 (g 2 + m 2 )) m| 

+ (g 4 + 5m 2 g 2 - Am% - (g 2 + m 2 ) m 2 ^) m| 

+ m^ (m^ - g 2 ) (2g 2 - 3m 2 + m 2 N ) m s + (g 2 - m 2 N ) 2 (m 2 N - m 2 ) , 
9o = m E (-2ra 2 m E g 4 - ( m N + m s) ( 3?T 4 _ 2 (Sm 2 ^ - 2m T m N + m|) m 2 

+ m^ (m^ - m 2 ) 2 (377% + m s )) g 2 + m N (m N - m s ) 2 m s (m^ + m s ) 3 ) 

= m N m^ - (g 2 + 2m 2 + m^) m| - m N (3g 2 - 3m 2 + 2m 2 N ) m| 
+ (g 4 - 5m 2 g 2 + 2m% + (g 2 + m 2 ) m^) m| 

+ m N (m 2 N - q 2 ) (2g 2 - 3m 2 + m 2 N ) m s + (g 2 - m 2 N ) 2 (m 2 - m 2 N ) , 
grW = (-m N q 6 + ((377% - m s ) m 2 + m^ {-2m 2 N + + m|)) g 4 

- m 2 (3m 2 - Am 2 N ) (m N - m s ) g 2 - m 2 (m^ - m s ) 3 (m N + m s ) 2 ) , 
/W = 3m N m% i + (g 2 - 2m 2 - Sm^) m^ + m N (3 (g 2 + m 2 ) - 4771^) m| 

- (g 4 + 5m 2 g 2 - 4777^ - (g 2 + m 2 ) m 2 N ) m| 

+ m N (m 2 N - q 2 ) (2q 2 - 3m 2 + m 2 N ) m 2 + (g 2 - m 2 N ) 2 (m 2 - m 2 N ) , 
g( b > = -m E (-2m 2 m s g 4 + (m N - m E ) (3m 4 - 2 (Sm^ + 2m E m 7V + m|) m 2 

+ m^v (3m N - m E ) (m^ + m s ) 2 ) g 2 + m^ (m^ - m s ) 3 m s (m^ + m s ) 2 ) 
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f- 



(«0 



f(c) 

Jo 



(c) 
3o 



/- 



(<9 



,(<*) 



-/o 



% 



m 2 (8m| + 5m N m% - (3g 2 + mfy) m| + 3m N [m 2 N - g 2 ) m s + (g 2 - m 2 N ) 2 



m 



N 



y-m N m% + (q 2 - 2m 2 N ) m| - Aq' 



m N m E 



- (g 4 + m 2 N q 2 - 2m%) m E + i 

- (m N - m E ) m E (m^ (2?^^ + m s ) g 4 + (m 4 - 2 (3m^ + 2m T m N + m|) m 2 

+ (mjy - ra s ) (m^ + m E ) 2 ) g 2 + 2m 2 (m^ + m E ) 2 (m 2 + m s (m E - m N ))) , 
m 2 ^8m E + 5m N m^ - (3g 2 + m 2 N ) 171^ + 3m N (m 2 N - q 2 ) m E + (g 2 - m%)' 

- (m N - m E ) 2 ^2m| - m N m% - (3g 2 + m 2 N ) m\ + 3™^ (m^ - g 2 ) m E 



(m N - m E ) m E ((m 4 - (2?^^ - Am j: m N - 2m|) m 2 + m^ (m N 



m 



N 



+ m E )) g 2 



+ (m N + m E ) 2 (2m 4 - 2 (2mfv - m T m N + m|) m 2 + m^ (77% - m E ) 2 (2772^ - m E ))) , 

(B5c) 



m 



iv 



m^ 



- (g 2 - 2m 2 N ) m| - 4g 2 m Ar m| 



+ (g 4 + m^g 2 - 2m^) m E + m^ (g 2 - m 2 N )' 

- m 2 (8m E - 5m N m^ - (3g 2 + m 2 N ) m| + 3m N (g 2 - mfy) m E + (g 2 - m 2 N ) 2 

m E (m^ + m E ) (-77% (277% - m E ) g 4 - (m 4 - 2 (3m 2 N - 2m Y m N + m|) m 2 

+ m^ (mjy - m E ) 2 (77% + m E )) g 2 - 2m 2 (m N - m E ) 2 (m 2 + m E (77% + m s ))) , 

(% + m E ) 2 (2m E + m^ml - (3g 2 + m%) m| + 3777^ (g 2 - m 2 N ) m E + (g 2 - m 2 N )' 



- m 2 (8ml - 5m N m% - (3g 2 + m^) m| + 3m N (g 2 - mfy) m E + (g 2 



m E (m^ + m E ) ((m 4 - 2 (m 2 N + 2m s m 7V - m|) m 2 + m^ (77% - m E ) (m N 
+ (77% - m E ) 2 (2m 4 - 2 (2m 2 v + m^rrij^ + m|) m 2 
+ m JV (m JV + m s ) 2 (2m JV + m s ))) . 



m s+ , m 



N 



\{m + m n ). 



m T 



m E ) J g 

(B5d) 
r0 ), the 



In our numerical computations, m E 
numbers being from Ref. [8]. 

In heavy baryon [14], the relevant Lagrangian can be found in Ref. [10], and the weak 

radiative and nonleptonic amplitudes in Eqs. (10) and (17) become, respectively, 



M(Bi — > Bp* 



-eG^B 



f 



- eG F B f 
M{Y> + -> Ntt) -- 



2(S-qS fl -S fl S-q)a + 2(S-qv fl -S fl v-q)b 



(q 2 - q^ v ■ q)c + 2 (g 2 - q" S ■ q)d 



B,,e: 



iG F mi + N A Nn + 2S ■ Vl 



B 



Ntt 



2m s 



(B6) 
(B7) 
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where v is the baryon four- velocity and S is the baryon spin operator. Following Ref. [10], to 
obtain the imaginary parts of the form factors we evaluate the loop diagrams displayed in Fig. 7. 
In the heavy-baryon approach, only the diagrams with the S + — > nir + transition yield nonzero 
contributions to the leading-order imaginary parts. The results are 



= P+ 5 )m - + ^^y/A^4[l + 



g 4 + 4m 2 (A 2 - g 2 ) 



4 (A 2 - g 2 ) 



3/2 



In 



A 2 - q 2 

2A 2 - q 2 - 2y/A 2 -ml^A 2 -q 2 

q A + 4 m 2^2 _ g2) 



(B8a) 



Im b 



(D + F) ml 
8V2 7rf n 



A. 



mt 



h 2 



A 3g 4 + 4m 2 (A 2 -g 2 ) ^ 



A 2 -g 2 
2A 2 



4(A 2 -g 2 ) 



5/2 



A 2 - 

2^/A ^ 



A 2 



g 4 + 4m 2 (A 2 - q 2 ) 



(B8b) 



Imc 



A (g 2 - 2m 2 ) 
+ g. i n 



A 2 - 2m 2 
A (A 2 - g 2 ) 



2(A 2 -q 2 ) 



3/2 



2A 2 - q 2 - 2 v /A 2 -m 2 v /A 2 -g 2 
g 4 + 4m 2 (A 2 - g 2 ) 



(B8c) 



T , (D + F)m 2 A 
Imd = - =^ — —A. 



a/A 



2 _ m 2 



§? 2 



* (A 2 -g 2 )' 



g 4 + 2g 2 A 2 + 4m 2 (A 2 -g 2 ) 
H ^ In 



4(A 2 -g 2 ) 



5/2 



2A 2 - g 2 - 2 v /A 2 -m 2 v /A 2 -g 2 
g 4 + 4m 2 (A 2 -g 2 ) 



(B8d) 



where A = m E —m N . We have checked that these formulas can be reproduced from the relativistic 
results in Eq. (B2) by expanding the latter in terms of A/m s , y / g 2 7m s , and m n /m j: and keeping 
the leading nonzero terms. 



s. y 



FIG. 7: Diagrams for imaginary part of S + — » £ry* amplitude. 
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APPENDIX C: REAL PARTS OF c(q 2 ) AND d{q 2 ) 

Vector mesons can contribute to c via the pole diagrams shown in Fig. 8(a). The strong vertices 
in the diagrams come from the Lagrangian [15, 16] 

c s = g D (BY{y,,B}) + g F (BY [\/ fl ,B}) + g (BYB)(y fl ) 

- \ e/ v ((2W - W) (?Q£ + (d IM A v - d w Aj , (CI) 

with V = |A 3 p° + • • • containing the nonet of vector-meson fields and D^\/ u = d^\J u + [V M , V y ] , 5 
whereas the weak vertices arise from 

£ w = G F m 2 n+ (h D (B{tfh£,B}) + h F (B [tfht,B]) + K(h^y^)) + H - c - , (C2) 

with h being a 3x3-matrix having elements h kl = 8 k2 8 zl which selects out s — > d transitions. The 
relevant parameters in C' s are Q D = —13.9 and Q F = 17.9 from a recent dispersive analysis [16, 17], 6 
and / v = 0.201 from p° -> e+e" rate [8], while those in £ w are h D = -72 MeV and h F = 179 MeV 
extracted at tree level from S-wave hyperon nonleptonic decays [18], but h v cannot be determined 
directly from data. To estimate hy, we use the SU(6)^ relation (7r°|7Y w |/i°) = (p°|7i w | A'* ) 
derived in Ref. [19]. Thus, employing the weak chiral Lagrangian = 7 8 f 2 (hd^Hd^^ + H.c, 



(a) 



(b) 

FIG. 8: Pole diagrams contributing to the c and d amplitudes. A single line (double line) denotes a baryon 
(vector meson) field, and a solid dot (hollow square) represents a strong (weak) vertex. 



5 Under a chiral transformation, V -► UVW and D^\l v -> UDWW. 

6 Although Q Q does not appear in our results, it enters the extraction of Q DF - Writing the pp\l part of C' s as 

hv^v {9 pN npI + g^NN^n + g^NN^n)' we have 9 p nn = D + Q F = 4 -°> 9 uN n = Q d + Q f + 2£ = 41 - 8 > and 

9<pNN = V 7 ^ (Q D — Qp + Q ) = —18.3, where the numbers are from Ref. [16, 17]. 
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with 7 8 = 7.8 x 10~ 8 from K — > tttt data, we find /i v = -4y 8 m 2 K / (G F m 2 n+ ) = -0.34 GeV 2 . 
Putting things together and adopting ideal oj-<j) mixing, we then obtain 

R ec = ^ v ( g ° ~ ^ m - + ^ D ~ kp } 



Q(m^-m N ) \q 2 -m 2 p q 2 - m* q 2 - m 2 

+ /y (go ~ gg) ^+ K ( 3 _ 1 + 2 \ 

12(g 2 - m 2 ^) yg 2 - m 2 g 2 - m 2 q 2 - m 2 J ' 

The form factor d can receive vector-meson contributions from the parity-violating Lagrangian 

C' w = G F rnZ r¥ h PV (ht{[B,'f 76 B],V ll }tf) + H.c. , (C4) 

which are represented by the diagram in Fig. 8(b). The parameter /i PV also cannot be fixed directly 
from data, and so we estimate it by adopting again the SU(6)^ results of Ref. [19] to be h PY = 2.41. 
It follows that 

Red =h^Zv(-A_--±_ + _l\ . (CB ) 
o \q mp q — q z — J 
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